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'■; ■ We have performed H and Ks band observations of the planetary system around HR 8799 using the new AO system 

(~ l ' at the Large Binocular Telescope and the PISCES Camera. The excellent instrument performance (Strehl ratios up 

PL(| to 80% in H band) enabled detection of the inner planet HR8799e in the H band for the first time. The H and Ks 

I , magnitudes of HR8799e are similar to those of planets c and d, with planet e slightly brighter. Therefore, HR8799e is 

P ■ likely slightly more massive than c and d. We also explored possible orbital configurations and their orbital stability. 

.4.^ ' We confirm that the orbits of planets b, c and e are consistent with being circular and coplanar; planet d should have 

t/^ I either an orbital eccentricity of about 0.1 or be non-coplanar with respect to 6 and c. Planet e can not be in circular 

, ^, , and coplanar orbit in a 4:2:1 mean motion resonances with c and d, while coplanar and circular orbits are allowed for a 
5:2 resonance. The analysis of dynamical stability shows that the system is highly unstable or chaotic when planetary 

fvj ' masses of about 5 Mj for b and 7 Mj for the other planets are adopted. Significant regions of dynamical stability for 

►^ I timescales of tens of Myr are found when adopting planetary masses of about 3.5, 5, 5, and 5 Mj for HR 8799 b, c, 

w-v , d, and e respectively. These masses are below the current estimates based on the stellar age (30 Myr) and theoretical 

f^ . models of substellar objects. 

r^ ■ Key words. (Stars:) individual: HR 8799 - Planetary systems - Instrumentation: adaptive optics - Techniques: high 

^^ ' angular resolution - Planets and satellites: dynamical evolution and stability - Planets and satellites: physical evolution 
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^ 1. Introduction planet (HR8799 e) at about 15 AU (jMarois et al.ll2010[ ). 
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The system is completed by three debris disk components. 

The planetary system around HR8799 represents a unique ^ ^^^^ ^f ^^^^ ^^^^ (y ^ ^50 j^) between about 6 to 10 

laboratory to constram the physical properties of massive ^U, a broad belt of cold dust (T ^ 45 K) between 90 

giant planets, to study the architecture of a crowded plane- ^^ ggg ^U, whose inner edge is probably defined by the 

H ■ tary system, and the hnk between planets and debris belts, interactions with the outer plan et, and a n extended halo 

^y: Three planetsjHR 8799 b, c and d) were discovered ^f s^j^^ll grains up to 1000 AU (ISu et al.l l200<^. The belt 

by| Maroisetal. | (120081), at a projected separation of about ^f ^^^^ just at about 1 00 AU ha ve been spatially resolved 

24, 38, and 68 AU, followed by the detection of an inner ^^ jq^^ ^^-^^ gpi^^^^. (|Su et al.ll2009l ). The central star is 

"" — ' — ~ ^ ^ an A 5 star located at 39.4 pc from the Sun (Ivan LeeuwenI 

Send offprint requests to: S. Esposito, [2^ characterized by A Boo-Uke abun dances anomalies 

e-mail: espositoSarcetri. astro, it i t^ 1 1 ,• rFi — t7 I N/^r>/^ K 

^ ^, ^\^^ . . , i. , „ , ,. • i-^ and 7 Doradus pulsations (IGrav fc Kavelll999l ). 

the LBl is an international collaboration among mstitu- ' ' 

tions in the United States, Italy and Germany. LBT Corporation 

partners are: The University of Arizona on behalf of the mi 1 -j. ^ r ^i TTno-TnA j. -ii --l r 

K . . .. , T i-i J- i\T • 1 J- A 4- c • Ihe architecture ot the HR8799 system, with its tour 

Arizona university system; fstituto Nazionale di Astrofasica, . , , , , , , , 1 , . r n 1 

Italy; LBT Beteiligungsgesellschaft, Germany, representing the S^ant planets and two belts, resembles that of our Solar sys- 

Max-Planck Society, the Astrophysical Institute Potsdam, and tem, especially when the two systems are plotted against 

Heidelberg University; The Ohio State University, and The the equilibrium temperature at various distances from 

Research Corporation, on behalf of The University of Notre the central star, taking the higher lu minosity of HR879 9 

Dame, University of Minnesota and University of Virginia. compared to the Sun into account (jMarois et al.l |2010[ ). 



S. Esposito et al.: LBT observations of the HR 8799 planetary system 



However, the planets around HR 8799 are much more mas- 
sive than those in the Solar System. 

The discovery of this planetary system prompted 
several investigations focused mos tly on the physi- 
cal properties of the planets (e.g. iBowler et al.l l2010t 
ICurrie et al.l 1201 it ) an d the dynamical stability of the 
planeta ry system (e. g. 'Gozd ziewski fc Migaszewskil l2009t 
rFabrvckv fc Murrav-Clav 20ld r 

The masses of the planets have not been determined 
dynamically and estimates are therefore deri ved from the 
stellar age of HR8799 and theoretical models. iMarois et al.l 
(|2008l) estimated age limits between 30 to 1 60 Myr, from 
the p osition of HR8799 on HR diagram. iMarois et al.l 
(|2010D narrowed the plausible age range to 30-60 Myr 
(with preference for the younger value) by classifying HR 
8799 as a probable member of the Columba moving group 
([Torres et all I2008D . The association of HR 8799 with 
Columba was questioned by iHinz et all (|2010D , who noted 
that the closest approach to the centroid of Columba mov- 
ing group was ~ 58 pc 27 Myr ago (this was considered 
too large for a direct link) and that the size of this and 
other young groups might be too large to have a common 
origin. Furthermore, available models of substellar objects 
suffer of significant uncertainties especially at young ages 
(jBaraffe et al.ll2002l ). leaving some ambiguity on the planet 
masses as derived from magnitudes or colors even at fixed 
age. 

On the other hand, studies of dynamical stability of the 
system showed that, for masses above 20 Mj, it is basically 
impossible to find orbital configurations compatible with 
the astrometric data and t hat are, at the same time, stable 
for the age of the system (|Moro-Martm et al.|[2010t) . The 
packed configuration of the system then favour the lowest 
planetary masses (corresponding to the youngest ages for 
the system), i.e. about 5 Mj for the outer planet and 7 Mj 
for the other three. Most of these studies were based on a 3- 
planet system architecture and should be extended includ- 
ing the fourth planet. The continuation of the astrometric 
monitoring and its "extension" to the past by identifica- 
tion o f some of the planets in improved reanalysis of past 
data (jSoummer et al.ll201lD is expected to provide tighter 
constrain on both the planetary orbits and masses. 

The development of new instrumentation for imaging of 
giant planets at small separation is crucial for a further un- 
derstanding of the system. Enhanced Strehl ratios allow to 
extend the detection space of the inner planets to additional 
wavelengths and then to better characterize their physical 
properties. In addition they allow to improve the accuracy 
of astrometric and photometric measurements thanks to 
the reduced speckle noise and the enhanced contrast of the 
planet PSF. 

The new adaptive optics system (FLAO) of the Large 
Binocular Telescope (LBT) has achieved, since its com- 
missioning, unprecedented p erformance with Stre hl Ratios 
higher than 80% in H band (JEsposito et al.l 120101) . Images 
of HR8799 system with such an instrumentation allowed 
us to detect for the first time the inner planet (HR8799e) 
in H band, enabling a more robust characterization of this 
planetary system. The present paper describes the instru- 
mentation and the observed procedures adopted to achieve 
this result, the data analysis procedures. Finally, it dis- 
cusses the results both in terms of physical characteriza- 
tion of the planets and of th e dynamical ar c hitect ure of the 
system. A companion paper. ISkemer et all (|2012[ ). presents 



3.3/im photometry of the HR 8799 system obtained using 
the Large Binocular Telescope Interferometer (LBTI) and 
the FLAO system, together with the independent analysis 
of the H band data discussed here and a comparison of the 
spectral energy distributions of the four planets with a va- 
riety of models including a new set of mixed-cloud models. 

2. Observations 

2.1. The LBT Adaptive Optics system 

The Large Binocular Telescope (LBT) is a unique telescope 
featurin g two co- m ounted optical trains with 8.4m primary 
mirrors (jHilll l2010l) . The First-Light Adaptive Optics sys- 
tem (FLAO) of the LBT takes advantage of two innova- 
tive key components, namely an adaptive secondary mirror 
with 672 actuators and a high-order pyramid wave-front 
sensor with a maximum pu pil sampling of 30 x 30 subaper- 
tures (jEsposito et al.ll2010l ). FLA0#1 system is located on 
the right telescope bent Gregorian focal station and con- 
trols the right secondary mirror. The WFS is mounted on 
rotating bearing that allows to compensate sky rotation: 
it receives the visible light reflected by a dichroic mirror, 
which transmits the infrared to the scientific instrument 
(in our case Pisces) . The commissioning of the instrument 
was completed in the winter of 2011, including a period of 
Science Demonstration Time (SDT) supplying a corrected 
wavefron t for the PISCES ima ging Near Infrared (NIR) 
Camera ([McCarthy et al.ll200lt) . The full adaptive optics 
imaging and spectroscopic channel will be completed when 
the LUCIFER (Mandel et al. 2000) instrument wfll be in- 
stalled on the NIR focal plane corrected by the FLA0#1. 
Since the initial phases of the on-sky commissioning, the 
FLA0#1 system reached performances never achieved be- 
fore on large ground-based optical telescopes. Images with 
40mas resolution and Strehl Ratios higher than 80% were 
obtained in H band (1.6 /zm). The images show a ratio be- 
tween the i ntensity at 0.4 arcsec and the central peak larger 
than lO"'' (JEsposito et al.ll201l[) . 



2.2. PISCES Camera 

The observations presented here used the 1-2.5 p,m cam- 
era PISCES (McCarthy et al. 200i|) at a bent- Gregorian 
focus of one 8.4 m primary mirror of the Large Binocular 
Telescope (LBT). Internal optics, cooled to 77 K, reimaged 
the f/15 focal plane onto a HAWAII-1 detector (1024 pixel 
square) at f/23.5, yielding a scale of 19.31 ± 0.03 mas/px 
with a field-of-view of 19.7 arcsec on a side (jClose et al.l 
l2012f ). This scale critically samples the diffraction-limit in 
the H-band (A/D = 40.5 mas). A cold pupil stop, nearly 
conjugated to the adaptive secondary mirror, shielded un- 
wanted background radiation. A dichroic beamsplitter lo- 
cated in the converging beam ahead of the camera directed 
visible light (< 0.95 p^m) onto the wavefront sensor unit. 
Images from PISCES are obtained from double-correlated 
sampling with a read-noise of 20 electrons. Further details 
on the performances of the PISCES camera coupl e d wit h 
the LBT AO system are described in iGuerra et al.l (|2012f ) . 



2.3. Observing strategy 

Our observations of HR 8799 were obtained with the Large 
Binocular Telescope in two different filters: H band on 
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Table 1. Summary of the observing setup and the observ- 
ing conditions for HR8799. 





H band filter 


Ks band filter 


Observation date (UT) 


201f Oct. f6 


2011 Nov. 09 


Number of images 


901 


328 


Exposure time 


2 s 


2 s 


Total integration time 


30 min 


11 min 


Total field rotation 


89.6° 


36.0° 


Seeing 


0.93 arcsec 


1.0 arcsec 



2011 October 16 UT, within the LBT PISCES+AO Science 
Verification Time (SV), and in Ks band on 2011 November 
09 UT during the Science Demonstration Time (SDT). The 
summary of the observing setup and of the observing con- 
ditions for these two epochs is reported in Table [T] 

The image rot ator was stopped to enable angular dif- 
ferential imaging (jMarois et al.l 120061 ). During the obser- 
vations in H band, in order to minimize the dynamic ef- 
fect of quasi-static speckles, we concentrated the observa- 
tions close to the meridian passage of the star achieving the 
largest angular coverage in the small time frame available. 
However, this strategy was not used during the observa- 
tion in Ks band because observations started slightly after 
the meridian passage due to technical problems. The seeing 
was measured with the LBT-DIMM pointing to the same 
direction of the scientific target on sky. 

We observed HR8799 saturating the inner region in both 
runs (i.e. in both filters H and Ks) at radii closer to 160 mas. 
At radii larger than 160 mas, the planets were observed 
with the stellar halo within the linear regime of the camera 
(that is below 8000-10000 counts). We collected calibration 
images to compute differential sky flat-field at the sunset 
and at the sunrise of the observing night. 



3. Data analysis procedure 

Initial image processing corrects all raw images for elec- 
tronic cross-talk between the quadrants in the detector us- 
ing Corquad, an IRAEQ task developed for this purposo 
The cross-talk coefficients were upd ated after final elec - 
tronic set-up at LBT as described in iGuerra et al.l (|2012[ ) 
After this, data analysis aiming to achieve the best possi- 
ble contrast was performed using two independent pipelines 
and data-reduction strategies called, hereafter " A" and " B" 
methods developed respectively at the Padova Observatory 
and at Steward Observatory. Both pipelines achieved co- 
herent results. 



3.1. Pipeline A 

The data analysis is composed of different steps 
that were implemented using IDL routines prepared 
for this purpose. The adopted procedure represents 
an optimization of data analysis routines prepared 
for similar datasets obtained with NACO at VLT 
(jChauvin fc The Naco Large Programme For Giant Planets 



^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 

■^ Corquad is available at jhttp : //ar ies . as . ar izona . edu/ 1 



l2010t ). The first step, which is critical because of heavy 
saturation, concerns the identification of the center of 
the star in each image. Accurate frame-to-frame relative 
positions were obtained by performing a cross-correlation 
(CC) between one image chosen as reference (the first 
image in the dataset) and all the other ones The peak of 
the CC represents the shift between the center position 
of the two images. The position of the maximum was 
obtained by means of a 2D Gaussian fit on this peak. This 
was found to be more precise than the Gaussian fit of the 
image. Much more difficult is to derive the position of 
the center of the star, with respect to that of faint planet 
images, that are not detectable in individual images. 
We proceeded as follows. First, the absolute position 
was obtained by finding the center of the stellar image 
after heavily smoothing using a 2D Gaussian fit. This 
procedure was repeated for every image. The results were 
then corrected for the relative frame-to-frame positions 
determined above, and then averaged. We found that on 
average the centers have an offset of about 0.5 pixels in 
both coordinates with respect to the value determined 
from the first image. We apply this correction to our star 
center position. The r.m.s. scatter for individual images 
is 0.31 pixels in both coordinates. Assuming that errors 
in the individual determinations are independent, a very 
small error of 0.01 pixels (~ 0.2 mas) is derived for the 
H-band data, and about twice this value for those in the 
K-band. However, systematic errors due to asymmetries 
of the PSF's are likely much larger. A rough estimate can 
be obtained by determining the star center in a slighly 
different way. Rather than heavily smoothing the profiles, 
we replaced pixel values in the saturated region with a 
constant value close to the maximum of the unsaturated 
pixels. The resulting center positions differ systematically 
by 0.116 ± 0.005 pixels with respect to the adopted ones. 
Since both procedure looks fairly legitimate, we conclude 
that our star position may have systematic errors as 
large as 0.1 pixels and likely more, as both procedures 
take into account only the outer regions of the PSF, and 
may underestimate asymmetries of these with respect 
to the core of the PSF, which is the only part that can 
be detected for the planets. Therefore, conservatively we 
adopt a larger error of 0.5 pixels in our discussion. 

In the second step, a 2D stellar profile is subtracted from 
each single image. The aim is to reduce the strong signal 
gradients present in the image, improving effectiveness of 
the following filtering procedures. This was done by sub- 
tracting to the values measured in each pixels the median 
of the counts in annuli one pixel wide at different separa- 
tions from the stellar center. After this, a low-pass filtering 
was applied to each image to eliminate bad and hot pixels. 
The next step was an high-pass filtering of the images. This 
was done by subtracting from each pixel the median of a 
sub-image composed by n x n pixels around it. The value of 
n changes according to the distance from the center of the 
star and it can be optimized in such a way to obtain the 
best final result. However, to reduce the possibility to self- 
-cubtvact j,i popciblc. companion, the median was performed 
-teSJSS tS'iS88Si^°flic central pixels of the sub-image. 

Finally, we implemented our v ersion of the angula r dif- 
ferential imaging (ADI) method (jMarois et al.ll2006D . We 
first located the center of the star of each image in the same 
position. Then, for each image, we selected a sub-sample of 
images in such a way that their rotation with respect to the 
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reference image was not too small, to avoid to self-subtract 
possible companion objects, and not too large, otherwise 
the speckle pattern would change too much and the im- 
age subtraction is no more effective in the speckle noise 
subtraction. The best criteria for image selection are dif- 
ferent at different separations from the central star, so we 
repeated this part of the procedure using different sets of 
images optimized for different separations. We then evalu- 
ated a median image for each sub-sample of images giving 
larger weights to images taken at shorter time lapses and 
the resulting image was then subtracted from the reference 
one. This subtracted image was then rotated by the appro- 
priate value, given by the parallactic angle as reported in 
the image header. To this angle, however, we had to add a 
further rotation to obtain a proper absolute orientation of 
each image. The calibration of this position angle is prelim- 
inary. However, an a posteriori test of the accuracy of the 
zero point correction to true north is given by the consis- 
tency of our astrometry with the predictions of the orbit of 
planet b by Soummer et al. (2011). This test indicates that 
the zero point is not in error by more than 0.3 degrees, in 
agreement with the error in the position of the true north 
by Close et al. (2012). This procedure is repeated for each 
image of the datacube. After this, the companion objects 
were in the same position in every image. We searched for 
such companions on the median of this datacube. 

For H band series, the best results in terms of the image 
S/N ratio for photometry at the observed angular separa- 
tion of HD8799e from the primary star are achieved us- 
ing only the frames taken closer to meridian passage (216 
frames corresponding to ±15°, for a total integration time 
of 7 minutes) . This is due to a better subtraction of speckle 
for larger angular velocity and shorter time baseline with re- 
spect to the full dataset. At larger separation, where speckle 
noise is less critical, the use of a larger dataset provides a 
slightly better S/N due to the larger photon flux. The re- 
sults presented in Sect. H] are based on the subset of the 216 
images taken close to meridian passage. 

For the case of the Ks band we used all the images 
(apart three of very poor quality) because it was not pos- 
sible to observe the star during meridian transit due to 
technical problems. 

3.2. Pipeline B 

The H-band data were independently reduced with a LOCT 
based pipeline (Locally Optimized Combination of Images; 
iLafreniere et al.ll2007b[) . which is hereafter referred to as 
Pipeline B. Pipeline B begins by flat/dark/distortion cor- 
recting the first 500 images (rather than 216 images for 
pipeline A), before the natural seeing worsened. Processed 
images were aligned by cross-correlation. The stellar pro- 
file was subtracted from each image, using an azimuthal 
average, and quasi-static speckles were suppressed by sub- 
tracting the median of the full set of images from each 
image. Through this step, pipeline A and pipeline B are 
mostly similar, other than the fact that pipeline B uses 
more data frames, and pipeline A includes an additional 
high-pass/low-pass filter (they are also independent im- 
plementations). After these (fairly standard) ADI steps, 
pipeline B uses LOCI to further suppress the noise of quasi- 
static speckle residuals. Pipeline B's results are statistically 
consistent with the results of Pipeline A. An in-depth de- 
scription of Pipeline B, and its photometric results are pre- 
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Fig. 1. Final image obtained from the H band data. 



sented in a companion paper, ISkemer et al.l (|2012[ ). In the 
context of this work, it is important to highlight that inde- 
pendent pipelines have produced similar results given that 
they are being used with a new instrumental setup. Some 
discrepancies are present in astrometric results between the 
two pipelines. This should probably arise from the different 
images centering procedure that, as said in the previous 
Paragraph, is probably the more tricky step of our data 
reduction procedures. 

3.3. Astrometric calibration and distortion correction 

A laser cut sieve mask was used in laboratory to derive the 
distortion correction for the PISCES camera (jGuerra et al.l 
.2012l ) . The coefficients of the polynomial computed in this 
way, reported in Appendix [^ have a 0.6- pixel accuracy in 
one a error. This error does not affect the data reduction 
since it produces negligible correction (less than 0.1 pix- 
els) on distances of the order of the separation of HR8799 
planets. 

The ab solute plat e sca le (19.274 mas/pix) was 
taken from IClose et al.l (J2012t) observations of the Orion 
Trapezium field. Preliminary analysis of on sky observa- 
tions of galactic Globular Cluster compared to HST cal- 
ibrated o nes further give s upport to the adopted value. 
Following IClose et all (|2012t) we also adopted a true north 
correction of 0.9 ± 0.3°. 



4. Results 

In Figure [1] and Figure [5] we display the final best ADI im- 
age obtained using the H- and K-band data, respectivelj0. 
In both figures the four planets are clearly visible; for best 
identification we marked their position with red circles. 
The FWHM of the planetary images are about 2.9 pixels, 
slightly larger than expected value of diffraction peaks. 

We measured both position and luminosity of the four 
detected planets. 

4.1. Astrometry 

In Table [5] and in Table |3] we report our astrometric mea- 
surements for the four planets around HR8799 in the H 



To reduce the noise visible in the final image, a smoothed im- 
age obtained using the IDL 'smooth' procedure was subtracted. 
Astrometry and photometry were performed on the original im- 
ages. 
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Fig. 2. Final image obtained from the K5 band data. 



Table 2. Astrometry (measured with respect to the cen- 
troid of the star HR8799) obtained from the H-band data 
(epoch 2011.79). 



Planet 


ARA (arcsec) 


ADec (arcsec) 


b 


1.579±0.011 


0.734±0.011 


c 


-0.561±0.010 


0.752 ±0.010 


d 


-0.299±0.010 


-0.563 ±0.010 


e 


-0.326±0.011 


-0.119 ±0.011 



Table 3. Astrometry (measured with respect to the cen- 
troid of the star HR8799) obtained from the Ks-band data 
(epoch 2011.86). 



Planet 


ARA (arcsec) 


ADec (arcsec) 


b 
c 
d 
e 


1.546 ±0.011 
-0.578 ±0.010 
-0.320±0.010 
-0.382±0.011 


0.725±0.011 
0.767±0.010 
-0.549±0.010 
-0.127±0.011 



Table 4. Astrometric error budgets (in mas) for the four 
HR8799 planets. 



Error source Error 



Notes 



and in the Ks-band respectively. Table S] gives the astro- 
metric error budget. We considered several sources of er- 
rors: centring of the star, image orientation (true north), 
scale distortion and planet centering errors (photometric 
errors). These last were evaluated comparing the measured 
position of a number of simulated planets, inserted into the 
original images, with their original values. Since these fake 
planets were inserted on the original images, these errors 
include artefacts due to data reduction and ADI, as well as 
the impact of speckles. 

Final errors in our astrometry were obtained by com- 
bining quadratically uncertainties concerning star position, 
true north, scale, distortion and statistical error in photom- 
etry. The dominant term is star centering. Uncertainties in 



Table 5. Photometry (absolute magnitudes) obtained from 
the H-band and the Kg-band data with the corresponding 
errors. Zero point photometric errors are not included. 



Planet 


H (mag) 


Ks (mag) 


b 


14.90 ±0.08 


13.98 ±0.06 


c 


13.90 ±0.12 


13.20 ±0.07 


d 


14.18 ±0.17 


13.71 ±0.25 


e 


13.53 ±0.43 


12.95 ±0.26 



Star pos. 


0.5 pix 


9 


9 


9 


9 


Same for all planets 


True north 


0.3 deg 


6 


4 


2 


1 


Error in pos. angle 


Scale 


0.2% 


2 


1 


1 


1 


Error in radial sep. 


Distortion 


0.15% 


2 


2 


2 


2 


Both in X and in y 


Phot. err. 










1 


5 


Both in X and in y 



the image orientation give a significant contribution for the 
most external planet while the photometric error is not neg- 
ligible for the planet e. 



4.2. Photometry 

Because of saturation of the center of the star and the lack 
of unsaturated exposures, we can not use the star magni- 
tude as a reference to evaluate the planet ones. For this 
reason we assumed, as a refe rence of our c alcula tion, the 
magnitude of HR87996 from iMarois et aO (|2008[ ) and we 
calculated the magnitude of the other planets by the ratio of 
planets counts. Pixels with separation < from the position 
of each planet were used for the relative photometry. We 
then repeated the same procedure assuming the same thing 
for HR8799c and we made a mean between the two results 
to obtain the final resuh. Given that HR8799& and HR8799c 
are located at quite large distances from the center of the 
star (1.72 and 0.97 arcsec respectively), this assumption is 
reasonable because the images of the planets are weakly af- 
fected by speckles. Table [5] reports the absolute magnitude 
of the four exo-planets as well as the associated errors. The 
magnitudes are corrected for the self-subtraction effect in- 
troduced by the ADI procedure. A set of errors (photomet- 
ric errors, error due to the self-subtraction of the SF of the 
star and errors due to the filtering data-reduction proce- 
dure) were evaluated by inserting into the original images 
~ 10 template planets with the same counts and at the 
same separation of each of the four planets. The same data 
analysis procedure was then performed on these latter im- 
ages. The standard deviation of the counts for planets at the 
same separation was then taken as the uncertainty on the 
photometry of each planet. Our H band pho tometry is on 
averag e about 0.20 mag fainter than that bv lSkemer et al.l 
(|2012 ). with a dispersion of 0.13 mag. This offset is primar- 
ily due to different ch oices for reference ma gnitude adopted 
for the outer p lanets, iMarois et al.l (l2008l) in o ur case and 
iMetchev etHu (|2009D for ISkemer et al.l (|2012D . which im- 
plies a systematic difference of 0.21 mag. 

Beside the agreeme nt with the photometry by 
IMarois et all (|2008i l2010l ) for planets b and c, which is ex- 
pected considering our normalization procedure, we note 
that the pho tometry of HR8799 e is also fully consistent 
with IMarois e t al. (2010) measurements. HR 8799(i results 
instead about 0.6 mag fainter in Ks band and 0.4 mag in H 
band. The discrepancy is however only marginally signifi- 
cant (about 2a) as HR8799(i happens to be projected close 
to the AO outer working angle in our images (see below). 
Therefore photometry of planet d is affected by larger errors 
(0.25 mag in Ks band). Further observations are necessary 
to confirm the reality and physical nature of this variation. 
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Fig. 3. 5ct contrast plot expressed in magnitude for the Fig. 4. 5cr contrast plot expressed in magnitude for the 

H-band data. Overplotted our own measurement for the Ks-band data. Overplotted our own measurement for the 

four HR8799 planet s (red triangle s) and the corresponding four HR8799 planet s (red trian gles) and the corresponding 

measurements from iMarois et al] (|2008.) as green squares. measurements from lMarois et al.. (2008. ) as green squares. 



4.3. Detection limits 

To calculate the detection limit as a function of the angular 
separation from the central star, we proceeded as follows. 
Along a radial direction, we calculated the standard de- 
viation of the intensity over a box oi d x d pixels (with 
d equal to the FWHM of the PSF i.e. 3 pixels for FLAP 
and PISCES) and one pixels step (jMasciadri et al.ll2005D . 
Values are then averaged over the azimuthal direction and 
the sum of the contrast AM at 5a versus the angular sep- 
aration is calculated normalizing with respect to the peak 
of the planet suitably re-scaled for the correct flux and tak- 
ing into account the flux losses due to the application of 
ADI. Fig. [3] and S] show the contrast obtained in H and 
Ks bands. The location of the four planets is marked with 
red triangles. For comparison , the position of the planets as 
determined bv lMarois et al.l (|2008i ) are showed using green 
squares. In Fig. [3] and S] the noise peak due to the outer 
working angle of the AO system (which is expected to be 
at 0.49 arcsec for H and 0.65 arcsec for kJ^ is clearly seen. 
Such a feature is not usually present in detection limits ob- 
tained with typical AO systems for 6-10 m class telescopes 
that achieve a much lower Strehl ratio than that delivered 
by the LBT FLAP system an d PISC ES camera with our 
observations (iMasciadri et al.ll2 005: La freniere et al.ll2007at 
iBiUer et al.ll2007l:lChauvhretal.ll2010D . 



'^ Following iRiccardi et al.l l|2010l ) the outer working angle 
(OWA) is computed as 

OWA = X/{2d) arcsec 



(1) 



where d is the effective inter-actuator distance considering the 
number of correcting modes and the actuator pitch projected to 
the primary mirror, which is 27 cm. Hence, d — 27 y^ 672 /n,nodes ■ 
When rimodes ~ 400 (which is the value in general used in 
FLAO), d = 35 cm. Hence we obtain OWA=0.49 arcsec (H 
band, 1.65 /im) and 0.65 arcsec (K band, 2.2 fim) It is impor- 
tant to mention that the extension of this halo spans about ± 
0.06 arcsec from th e OWA values reported above (Fig. 17 in 
IRiccardi et al.lf2010l ). These values agree very well with the po- 
sition of the secondary maxima in the contrast vs separation 
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Fig. 5. H/H— K color magnitude diagram comparing the 
observed magnitudes of IIR8799 planets with other cool 
objects in the field and known low mass companion. Three 
planets of the HR8799 system (c,d,e) are quite consistent 
with the observed L and T spectral type BDs. IIR8799 b is 
manifestly redder and fainter than L and T spectral type 
BDs separation. 



5. Physical properties of planets around HR8799 

We compared the near-IR properties of all HR8799 plan- 
ets with those of other substellar structures in order to 
investigate systematic differences due to e.g. mass and age. 
First of all we compared the H and Ks magni tude observed 
for IIR 8799 planets with those obtained by iLeggett et al.l 
(|2010D for field brown dwarfs (BDs) . We then collected from 
the literature data for others low mass companions with di- 
rect imaging observations in the same photometric bands 
(Table IB]). All these objects, together with the planets of 
HR8799, may be plotted in color-magnitude diagrams built 
with H and Ks filters. 

The color-magnitude diagrams is shown in Figure O 
Three planets of HR8799 (cde) have colors roughly con- 
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Table 6. Photometry of others smaU mass companions adopted by hterature. All listed magnitude are absolute magni- 
tude. 



Companion 



D 

pc 



Mass 
Mj 



Age 
Myr 



H 



K 



Ref. 



2M1207b 


52.4 


5 


8 


16.40 ± 0.20 


14.49 ± 0.21 


13.33 ± 0.11 


1,2 


lRXJ1609.1-210524b 


140.0 


8 


5 


12.17 ± 0.12 


11.14 ± 0.07 


10.44 ± 0.18 


3,4 


AB Picb 


47.3 


13 


30 


12.80 ± 0.10 


11.31 ± 0.10 


10.76 ± 0.08 


5 


HD203030b 


40.8 


23 


130-400 


15.08 ± 0.55 


13.80 ± 0.12 


13.16 ± 0.10 


6 


HIP 78530b 


156.7 


23 


5 


8.94 ± 0.23 


8.33 ± 0.22 


8.14 ± 0.22 


7 


CD 35 2722b 


21.0 


31 


100 


11.99 ± 0.18 


11.14 ± 0.19 


10.37 ± 0.16 


8 


SR 12 AB c 


125.0 






9.9 ± 0.7 


8.8 ± 0.9 


8.6 ± 1.1 


9 


HN Peg b 


18.4 


21 


300 


14.54 ± 0.06 


14.08 ± 0.06 


13.75± 0.06 


10 


Ross458 (AB) c 


11.7 


8.5 


475 


16.42± 0.67 


16.78 ± 0.69 


16.50 ± 0.69 


11 


CSC 06214 -00210 b 


145.0 


17 


11 


10.5 


9.6 


9.1 


12 


2M 044144 b 


140.0 


7.5 


1 




9.89±0.10 


9.21 ± 0.10 


13 


CFBDSIR J1458-H013b 


23.1 


6.5 


3000 


19.84 ±0.40 


20.69 ± 0.27 


21.02 ± 0.37 


14 


DH Tau b 


143.5 


11 


1 


9.65±0.05 


9.03±0.04 


8.33±0.02 


15 



Refer ences. References: fP): 'Chauvin et al.l ||2004 D: (2) : 'Mohanty ct al.' ('2007') (3): 'Lafreniere et al' (^ 

Il201iri: ('5'):[C hauvin et al. (2005); (6): M etchev fc H illenbrand (2006); (7): Lafreniere ct al. (2011); (8) _^ ^ 

iKuzuhara eTa l. (2011) (lOlilLuhman et aP f2007): (ll): IBurningliam et al.l (|2011l ): (12): [Ireland et all (120111) : (13p lTodorov et all 
l|2010H : fl4)T ]Liu et aP (J201lh : fl5'): l"toli et al., (,2005h 



: (A): [Ireland et all 
Wa hhai et al.1 (1201111: (9) 



sistent with L and T BDs sequence. In particular HR8799c 
and HR8799(i are slightly underluminous with reference to 
the rim between L and T boundary, with HR8799(i red- 
der than the colors of T BDs. We also recall that our Ks 
phot ometry of d has rat her large error and is 0.6 fainter 
than [Marois et all ()2008[ ) one (see Sect. 15). iMarois et al] 



6. The architecture of the HR8799 planetary 
system 

6.1. Orbital fit 



(|2008D photometry would imply a redder color, and a po 
sition slightly outside the field BD sequence. HR87996 has 
redder colors with respect the others HR8799 planets, and 
it is manifestly under luminous with respect not only the 
L-T sequence of BDs but also with respect to the other 
planets of the HR8799 system (regardless of the systematic 
uncertainties in the photometry). A similar but even more 
extreme anomaly is that of 2M1207b. To explain these lo- 
cations in the color-magnitude phase space various ad-hoc 
hypothesis were proposed, such as an occulting edge-on cir- 
cumplanetary disk (Mohanty et al. 2007) or a collision af- 
terglow (JMamaiek fc Meveril2007^ . The sim ilar discrepan- 
cies oc curring for HR87996 and 2M12 07b (iBarman et~al] 
I2OIIHI and specific inconsistencies (see ISkemer et al.ll201ll 
for details) argue against such ad-hoc explanations suggest- 
ing a more general feature linked to cloud properties in 
low-gravity atmospheres. 



In summary, we can confirm the faintness and redder 
NIR colors of HR8799b. When considering only H/H-K 
color-magnitude diagram, the positions of HR8799c,(i, and 
e is close to the field brown dwarf sequence, especially for 
the brightest planet e. Including information from other 
wavelengths clearly suggests the HR8799c, d, and c show 
significant differences. This is likely due to the presence 
of thick cloud laye rs and no n equilibrium chemi stry in 
their atmospheres (Currie et al. 201 P Barman ct al. i i2011at 
ISkemer et al.ll2013) . A more extensive disc ussion is given in 
the companion paper (jSkemer et al.l[2012n . 



Since the planet discovery papers (jMarois et al.l I2008L 
I2O1O0 . it resulted that all the planets in the HR8799 sys- 
tem orbit the star in the same direction and their orbits are 
roughly compatible with nearly pole-on circ ular orbits. The 
cumu lation of additional observational data (jBergfors et al.l 
l2011i ) and the detection of the some of pla nets in older im- 
ages, thanks to a repr ocessi ng of the data (ISoumrner et al.l 
I201li iFukagawa et al.ll2009l: iLafreniere et al.ll2009D allowed 
to further extend the time baseline of the observations. 

The time coverage from 1998 to present epoch rep- 
resents only a minor fraction of the orbital periods. 
Nevertheless some indications of the actual orbits is emerg- 
ing from the accumulated data. The orbits of planets b and 
c are compatible with circular orbits seen nearly but not 
exactly pole-on. Instead the orbit of planet d is eccentri c 
and/or seen at a different inclination ("Bcrgf ors et al.l[201ll ). 

A recent study of the configuration o f the outer part 
of the HR 8799 system was presented by ISoummer et al.l 
(|2011l ). taking advantage of the detection of the three outer 
planets in HST images taken in 1998. They consider as the 
most likely solution a coplanar system in 4:2:1 mean motion 
resonance, with HR8799d in a slightly eccentric orbit (e — 
0.1), as suggested in the dynam ical stability analysis by 
iFabrvckv fc Murrav-Cla^ (1201011 . 

The increase of the time baseline of 1 yr is not de- 
cisive for a major revision of these results. However, 
ISoummer et al.l (|2011l ) consider only the outer three planets 
and they restricted to coplanar orbits. Furthermore, evalu- 
ation of dynamical stability of the proposed configuration 
was not performed. 

Fig. in] shows th e relative astrometry of HR 8799 plan- 
ets, as compiled by Bcrgfors et al. (201Y) pl us the in c lusion 
of data from Galicher et al. (2011), Soum mer et al.l (J201l|) 
and our own measti r ement s (Tables Hj and |3|) . Overplotted 
the ISoummer et al.l (|201ll) orbital solutions for the three 
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Fig. 6. Relative astrometry of HR 8799 planets. Overline are the ISoummer et al.l (|2011[ ) orbital solution for the outer 
planets and a coplanar circular orbital solution in 5:2 resonance for the inner planet (blue lines, Case A in Tabic [7]). For 
HR8799e, the dotted line shows the circular and coplanar orbital solution for 2:1 mean motion resonance with planet 
d, which does not fit the observations. Filled circles: literature results. Filled squares: data from the present paper. Red 
lines connect the predicted and observed position for all the data points. 
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Fig. 7. Residuals in x and y coordinates vs time for each of the planets with respect to the orbital solution plotted in 
Fig.EKCase AinTablelZl). 



outer planets and our own solution for the inner planet 
(see below for the fitting procedure). 

New data are consistent with t he previous works and 
follow the best-fit orbit derived in ISoumnier et al] (12011 ) 
for the three outer planets. If we adopt the lSoummer et all 
(|2011f) orbital solution, planet e can not be in a circular and 
coplanar orbit with a 2:1 mean motion resonance with d. 
Instead, a 5:2 resonance represents a satisfactory fit to the 
data, when assuming circular orbit and coplanarity with 
the other planets (Fig. El Case A in Table [T]) . 

This orbital solution represents only one of the possi- 
ble ones fitting the observational data. Unique determina- 
tion of orbital parameters is not yet possible as the ob- 
servations cover only a minor fract ion of the orbital peri - 
ods. To somewhat complement the ISoummer et al.l (|2011[ ) 
study, that assumes coplanarity between the three outer 
planets, we focused our attention here on non-coplanar con- 
figurations. We restrict our analysis to circular orbits to 
reduce the number of parameters. However, this approxi- 
mation is likely not realistic considering that secular per- 
turbations among planets causes some eccentricity pump- 
ing, with average values of about 0.03-0.05. We also con- 
si dered only a stellar ma ss of 1.56 Mq, the preferred value 
in ISoummer et al.l ()2011| ) . Our analysis of possible orbital 
solutions was performed in two steps: 



— MonteCarlo simulations for a broad exploration of the 
orbital solutions compatible with the data and corre- 



lations between the orbital parameters (Fig. |5|). This 
allows us to identify the possible ranges to bound the 
least-square orbital fitting and to identify appropriate 
initial guesses. 
— Best-fit least-square orbital solution for the four plan- 
ets simultaneously using the Levenberg-Marquard min- 
imization algorithm as implemented in the IDL routine 
MPFIT. The program allows to fix some of the param- 
eters (e.g. impose null eccentricity) or to tie some pa- 
rameters of one planet to those of another planet (e.g 
for coplanar orbits or for imposing orbital resonances). 

When forcing all consecutive planets to be involved in 
2:1 resonance (8:4:2:1 resonance involving all the known 
planets) we found orbital solutions that nicely fit the avail- 
able data. The formal best-fit case is listed as case B in 
Table [T] A notable property of this orbit is the strong de- 
viation from coplanarity involving the inner planet e, larger 
than 50° with respect to planet d (when including both i 
and fi into account). Forcing planets bed to be in a 4:2:1 
resonance, planets de to be in a 5:2 resonance, and assum- 
ing all orbits are circular yields the orbital parameters la- 
beled as case C in Table [71 In this case the deviations from 
coplanarity are or the order of 10 — 20° and can become 
much smaller when allowing for some eccentricity of the 
orbits (see Case A). Such differences in the orientation of 
HR8799e orbit depending on the chosen resonance with d 
are independent on the orbit of planet b. We found a similar 
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Fig. 8. Reduced x^ surfaces for orbital period, inclination and longitude of node for HR8799 planets assuming circular 
orbits. Each plot row refers to one of the planets, starting from b (top) to e (down). The plotted levels refer to reduced 
X values in step of 0.20 starting from 2.0, with darker area correspinding to lower x values. 



behaviour when adopting an outer orbit for b, oii t side 2 :1 
resonance with c following ISudol fc Haghighipoud (|2012t ). 

Additional constraints on the orientation of the compo- 
nents of the H R 8799 system are a vailable for the central 
star, for which I Wright et all (|201ltl derived i > 40° from 
asteroscismology, and for the outer debris disk, for which 
ISu e t al. (2009) found i < 25°. The best-fit solutions listed 
in Table [7] show that the outer three planets have inter- 
mediate values, with hints of a trend in inclination from b 
to d toward more pole-on orientations inside (cases B and 
C). The relative inclination of e instead depends critically 
on the adopted orbital period, as discussed above. The two 
orbital solutions (B and C) while fitting in the same way 
the current data will diverge significantly within 2-3 years. 



6.2. Dynamical stability 

After these explorations based only on observational data, 
we consider the additional clues that can be derived from 
the dynamical stability of the system. 

Previous dynamical studies of the system considering 
the three planet systems (without planet e), found that only 
a minor fraction of the orbital solutions compatible with the 
astrometric observations is dynamically stable. Stability is 
favoured by the occurrence of a 2:1 mean motion orbital 
resonance between each couple of consecutive planets, low 
eccentricity of their orbits and by planetary masses at the 
low end of the plausible values derived from their lumi- 
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Fig. 9. Relative astrometry of HR 8799 planets available in literature. Blue lines show the adopted orbital solution (Case 
B in Table [71 non coplanar circular orbits). Red lines connect the predicted and observed position for all the data points. 
For HR8799e, the dotted line shows the orbital solution for 5:2 mean motion resonance with planet d (Case C in Table 
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Fig. 10. Residuals in x and y coordinates vs time for each of the planets with respect to the orbital solution plotted in 
Fig.|9](Case Bin Table [T]). 



nosity and stellar age (iQozdzie wski &: Migaszewski 2009t 
iReidemeister et al.| [2009: Fabryck 7fc'Murrav-Clavll2010[ ). 



Recently, ISudol fc Haehighipourj (|2012[ ) presented a 
study of the dynamical stability of the 4-planet system. 
They found that the system is strongly unstable with a few 
solutions that barely survive for ages comparable to the 
age of the syst em. Planet e is marginally fit by these stable 
solutions, and lSudol fc Haghighipoun (|2012l) proposed that 
HR8799e astrometry is systematically biased toward larger 
separations. Additional, more limited, dynamical studies of 
the 4- planet systern were perforr ned in the discover y pa- 
per bv lMaroiset~an (|2010t) and in lCurrie etaP (|201lt) with 
similar results. 

We numerically integrated the orbits of the planets 
around the star in a full 5 body problem using the nu- 
merical integrator RADAU (JEverh art 1985). The timespan 
of the integration was fixed to 25 Myr and the computed 
orbital elements were used to estimate the stability propcr- 
ti es using the Freauancy Map Analysis (FMA) technique as 
in iMarzari et al.l (200o). We studied the shift in frequency 
of the perihelion longitude of the middle planet as stability 
indicator. The phase space of a four planet system is wide 
so we devised the following strategy. We set the two middle 
planets in a 2:1 resonance while the outer planet is moved 
from the 2:1 to the 7:3 resonances with the third planet. 
Planet e is started with random orbital elements close to 
its expected orbit for a 2:1 and 5:2 mean motion resonances 
with planet d. The orbital elements are sampled randomly 



around this initial configuration looking for solutions which 
may be stable for long timespans. 

Assuming a mas s for the star of 1.56 Mq 
(jSoummer et al.l 1201 ll ) and for the planets the nomi- 
nal masses 9,10,10,7 Mj for e,d,c,b, respectively and by 
numerically integrating the orbits we find highly chaotic 
solutions that survive on orbits similar to the observed 
ones for a few Myrs only. This is not compatible with the 
stellar age even assuming that the primordial gas disk 
damped the eccentricities of the planets granting stability. 
The disk lifetime can be no longer than ^^ 20 x 10^ 
years. This is a conservative assumption as, according to 
observations, the presence of gas in the inner regions of the 
disk is observed o nly for timescales shorte r than 10^ years 
(see the review bv lWilliams fc Ciezall2011l) . Of course, due 
to the chaotic nature of the orbits, it may be possible to 
find orb its surviving for longer t imescales, as those found 
by .Sudol fc Haghighipouij ()2012t ). but a different choice of 
the numerical integrator gives already different solutions 
undergoing close encounters a fter a f ew M yrs. As already 
noted by Fabrycky & Murra v-Clavl (|2010f ). reducing the 
planetary masses leads to longer timespans before the 
onset of a 'Jumping Jupiter' phase. By reducing all masses 
by 30% and setting the outer planet in a 7:3 resonance 
with the third planet and planet e in a 5:2 resonance with 
the second planet, we find the most stable solutions to 
the system, according to the FMA analysis. They are still 
chaotic but they survive at least over the estimated stellar 
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Table 7. Orbital solu tion considered in this p aper. Case A: 
Orbital solution from ISoummer e t al.l (|201l[ ) + our fit for 
planet e (circular orbit, coplanar with other planets and in 
5:2 mean motion resonance with planet d). Case B: best fit 
orbital solution imposing circular orbits and 8:4:2:1 mean 
motion resonance. Case C: best fit orbital solution imposing 
circular orbits and 4:2:1 mean motion resonance for planets 
bed and 5:2 for de 



are 2, 6, and 3 Mj and ages of 10, 20 and 10 Myr for planets 
6, c, and d, respectively. 

Our explorative analysis can not be considered com- 
pleted. For example, we did not explored the possible im- 
pact of varying the stellar mas s within the ra t her br oad er- 
rorbars (±0.3 Mq) allowed by iGrav fc Kavd ()1999[ ) study. 
A full analysis is postponed to a forthcoming paper. 



Parameter 


A 


B 


C 


A (yr) 

4(°) 

n, (°) 


449.7 
28.0 
35.5 


448.44 
20.07 
22.52 


448.44 
20.07 
22.53 


et 

TOb (yr) 
at (AU) 


1997.55 
68.08 


2015.37 
67.94 


2015.37 
67.94 


P. (yr) 

ic (deg) 


224.9 
28.0 
35.5 


224.21 
27.78 
45.86 


224.22 
27.78 
45.86 


Gc 
^c(°) 

TOc (yr) 
a, (AU) 


1844.10 
42.89 


1838.89 
42.80 


1838.89 
42.80 


Pd (yr) 
idC) 

fid 

^d(°) 

rOd (yr) 

ffld (AU) 


112.4 
28.0 
35.5 
0.10 
80.2 
1992.31 

27.01 


112.11 
40.43 
36.76 

1965.20 
26.97 


112.11 
40.43 
36.76 

1965.20 
26.97 


Pe (yr) 


44.96 
28.0 
35.5 


56.05 
31.13 
-62.81 


44.84 
20.09 
43.39 


fie 
^e(°) 

TOe (yr) 

fle (AU) 


1987.08 
14.66 


1996.70 
16.99 


1986.14 
14.64 


Mstar (Mq) 


1.56 


1.56 


1.56 



age. However, these solutions are not fully compatible 
with the observational data. If the mass of the planets is 
reduced at 50% of the original values, the phase space area 
where orbits survive for at least 30 Myr becomes significant 
and even solutions with all four planets in mutual 2:1 
resonance become stable over a longer timspan. 

While a 30% reduction of the planetary masses is still 
compatible with the stellar age and theoretical models, a 
reduction of 50% (about 3.5, 5, 5, and 5 Afj for HR 8799 
b,c,d,e respec ti vely) is below the current estimates. Using 
iBaraffe et aD (|2003[ ) models, a stellar age of 15-20 Myr 
would be required to be compatible with such low masses 
for the plane t , whic h is below the youngest age derived by 
iMarois et al.l (J2010t) from membership to Columba associa- 
tion (nominal age 30 Myr) and only marginally compatible 
with t he younges t age estima t e (20 Myr) by I Moor et al.l 
(|2006D . However. [Torres et al.l (|2008[ ) quoted a significant 
age uncertainty for Columba association. Furthermore, 
some age dispersion within the association or a small age 
difference between the star and the planets might help to 
reconciliate evolutionary model predictions and constraints 
from dynamical stability of the system. The 3a lower limits 
of the model atmospheres fit bv lMadhusudhan et al.l (|201l[ ) 



7. Conclusion 

We have performed H and Ks band observations using the 
new AO system at the Large Binocular Telescope and the 
PISCES Camera. Analyses performed independently by dif- 
ferent members of our team using two different pipelines 
yielded consistent results. The excellent instrument perfor- 
mance (Strehl ratios up to 80% in H band) enabled de- 
tection of the inner planet HR8799e for the first time in 
the H band. The H and Ks magnitudes of HR8799e are 
similar to those of planets c and d, with planet e slightly 
brighter. Therefore, the inner planet should have similar 
masses to c and d, being likely slightly more massive. When 
considering only H/H— K and K/H— K color-magnitude dia- 
grams, the positions of HR8799c,d, and e is close to the field 
brown dwarf sequence, especially for the brightest planet 
e. Planet b is instead significantly fainter than the other 
HR8799 planets, as already discussed in the literature. The 
occurrence of some systematic differences at 0.2 mag level 
between photometric measurement of the HR8799 planets 
from different sources have some impact on these conclu- 
sions and should be further investigated. 

We also collected the available astrometric measure- 
ments of HR8799 planets, exploring possible orbital con- 
figurations and their orbital stability. We confirm that the 
orbits of planets 6, c and e are consistent with being cir- 
cular and coplanar; planet d should have either an orbital 
eccentricity of about 0.1 or be non-coplanar with respect to 
b and c. We found that the planet e can not be in a circular 
and coplanar orbit with the other planets and in 4:2:1 mean 
motion resonances with planets c and d. These resonances 
require significant deviations from coplanarity or eccentric- 
ity. A coplanar and circular orbit with 5:2 resonance be- 
tween d and e is instead compatible with the observational 
data. 

We found the system to be highly unstable or chaotic 
when the nominal planetary masses are adopted. Significant 
regions of dynamical stability for timescales of tens of Myr 
are found only when adopting planetary masses of about 
3.5, 5, 5, and 5 Mj for HR 8799 6, c, d, and e respec- 
tively. These masses are below the current estimates based 
on stellar age (30 Myr) and theoretical models of substel- 
lar objects. A more complete exploration of the parameters 
space will be performed in a forthcoming study to iden- 
tify possible system configurations that are compatible with 
the observations and dynamically stable. In any case it is 
more difficult to find dynamically stable solutions for the 4- 
planet system than for the 3-planet case considered in most 
of the literature studies. This implies smaller upper limits 
on planetary masses from dynamical stability constraints. 
These results will be relevant to place clues on the physical 
models of planet structure and atmospheres. 

On the observational side, the continuation of the astro- 
metric monitoring is mandatory for a better characteriza- 
tion of the system. With an orbital period of about 50 yr for 
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HR8799e, a few more years of observations will allow signif- 
icant refinement of the constraints we can put on the orbital 
parameters. New observations should pay specific attention 
to achieve the best astrometric accuracy (optimization of 
the observing procedure and instrument set-up, dedicated 
astrometric calibrations, etc.). 

Appendix A: Distortion correction coefficients 

Drizzle coefficients for PISCES distortion and plate scale 
determination were obtained for drizzled and un-drizzled 
PISCES images from the sieve mask data. For what con- 
cerns the fitted Drizzle coefficients, let us assume x' ,y' are 
corrected centroid values in pixels, x, y are raw data cen- 
troid values, in pixels. Then xO, yO translate the distortion 
equation to an appropriate centre for the distortion equa- 
tion: 

x' = aa + ai{x — xQ) + a^ix — xQ)^ + a(i{x — a;0)'^+ 

a2(z/ - yO) + a^{x - x{)){y - yQ) + a'j{x - x{))\y - yO)4 
a5(?/ - yO)2 + as{x - xO){y ~ yOf + a^iy - yOf 

y' = bo + bi{x - xO) + 63(0; - xOf + be{x - xOf+ 

hiy - yO) + bzix - xO){y - yO) + br{x - xOf{y - yO)+ 
b^{y - yOf + bs{x- xQ){y - yOf + 69(2/ - yOf 

The values of the ai and hi coefficients are listed in the 
Table EH 
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